Search for: All records

The parameterization of subgrid‐scale processes such as boundary layer (PBL) turbulence introduces uncertainty in Earth System Model (ESM) results. This uncertainty can contribute to or exacerbate existing biases in representing key physical processes. This study analyzes the influence of tunable parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBBX) scheme. CLUBB is the operational PBL parameterization in the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community ESM version 2 (CESM2). We perform the Morris one‐at‐a‐time (MOAT) parameter sensitivity analysis using short‐term (3‐day), initialized hindcasts of CAM6‐CLUBBX with 24 unique initial conditions. Several input parameters modulating vertical momentum flux appear most influential for various regionally‐averaged quantities, namely surface stress and shortwave cloud forcing (SWCF). These parameter sensitivities have a spatial dependence, with parameters governing momentum flux most influential in regions of high vertical wind shear (e.g., the mid‐latitude storm tracks). We next evaluate several experimental 20‐year simulations of CAM6‐CLUBBX with targeted parameter perturbations. We find that parameter perturbations produce similar physical mechanisms in both short‐term and long‐term simulations, but these physical responses can be muted due to nonlinear feedbacks manifesting over time scales longer than 3 days, thus causing differences in how output metrics respond in the long‐term simulations. Analysis of turbulent fluxes in CLUBBX indicates that the influential parameters affect vertical fluxes of heat, moisture, and momentum, providing physical pathways for the sensitivities identified in this study. 

Improving the prediction of clouds in shallow-cumulus regimes via turbulence parameterization in the planetary boundary layer (PBL) will likely increase the global skill of global climate models (GCMs) because this cloud regime is common over tropical oceans where low-cloud fraction has a large impact on Earth's radiative budget. This study attempts to improve the prediction of PBL structure in tropical trade wind regimes in the Community Atmosphere Model (CAM) by updating its formulation of momentum flux in CLUBB (Cloud Layers Unified by Binormals), which currently does not by default allow for upgradient momentum fluxes. Hindcast CAM output from custom CLUBB configurations which permit countergradient momentum fluxes are compared to in situ observations from weather balloons collected during the ElUcidating the RolE of Cloud–Circulation Coupling in ClimAte and Atlantic Tradewind Ocean–Atmosphere Mesoscale Interaction Campaign (EUREC4A/ATOMIC) field campaign in the tropical Atlantic in early 2020. Comparing a version with CAM–CLUBB with a prognostic treatment of momentum fluxes results in vertical profiles that better match large-eddy simulation results. Countergradient fluxes are frequently simulated between 950 and 850 hPa over the EUREC4A/ATOMIC period in CAM–CLUBB. Further modification to the planetary boundary layer (PBL) parameterization by implementing a more generalized calculation of the turbulent length scale reduces model bias and root mean squared error (RMSE) relative to sounding data when coupled with the prognostic momentum configuration. Benefits are also seen in the diurnal cycle, although more systematic model errors persist. A cursory budget analysis suggests the buoyant production of momentum fluxes, both above and below the jet maximum, significantly contributes to the frequency and depth of countergradient vertical momentum fluxes in the study region. This paper provides evidence that higher-order turbulence parameterizations may offer pathways for improving the simulation of trade wind regimes in global models, particularly when evaluated in a process study framework. 

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models. 

Recent studies have demonstrated that high-resolution (∼25 km) Earth System Models (ESMs) have the potential to skillfully predict tropical cyclone (TC) occurrence and intensity. However, biases in ESM TCs still exist, largely due to the need to parameterize processes such as boundary layer (PBL) turbulence. Building on past studies, we hypothesize that the depiction of the TC PBL in ESMs is sensitive to the configuration of the PBL parameterization scheme, and that the targeted perturbation of tunable parameters can reduce biases. The Morris one-at-a-time (MOAT) method is implemented to assess the sensitivity of the TC PBL to tunable parameters in the PBL scheme in an idealized configuration of the Community Atmosphere Model, version 6 (CAM6). The MOAT method objectively identifies several parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBB) scheme that appreciably influence the structure of the TC PBL. We then perturb the parameters identified by the MOAT method within a suite of CAM6 ensemble simulations and find a reduction in model biases compared to observations and a high-resolution, cloud-resolving model. We demonstrate that the high-sensitivity parameters are tied to PBL processes that reduce turbulent mixing and effective eddy diffusivity, and that in CAM6 these parameters alter the TC PBL in a manner consistent with past modeling studies. In this way, we provide an initial identification of process-based input parameters that, when altered, have the potential to improve TC predictions by ESMs. 

Abstract. This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs. 

Abstract Nudging is a ubiquitous capability of numerical weather and climate models that is widely used in a variety of applications (e.g., crude data assimilation, “intelligent” interpolation between analysis times, constraining flow in tracer advection/diffusion simulations). Here, the focus is on the momentum nudging tendencies themselves, rather than the atmospheric state that results from application of the method. The initial intent was to interpret these tendencies as a quantitative estimate of model error (net parameterization error in particular). However, it was found that nudging tendencies depend strongly on the nudging time scale chosen, which is the primary result presented here. Reducing the nudging time scale reduces the difference between the model state and the target state, but much less so than the reduction in the nudging time scale, resulting in increased nudging tendencies. The dynamical core, in particular, appears to increasingly oppose nudging tendencies as the nudging time scale is reduced. A heuristic analysis suggests such a result should be expected as long as the state the model is trying to achieve differs from the target state, regardless of the type of target state (e.g., a reanalysis, another model). These results suggest nudging tendencies cannot bequantitativelyinterpreted as model error. Still, two experiments aimed at seeing how nudging can identify a withheld parameterization suggest nudging tendencies do contain some information on model errors and/or missing physical processes and still might be useful in model development and tuning, even if only qualitatively.